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We study the Fermi energy level dependence of nonlinear terahertz (THz) transmission of gated 
multi-layer and single-layer graphene transferred onto sapphire and quartz substrates. The two 
samples represent two limits of low-field impurity scattering: short-range neutral and long-range 
charged impurity scattering, respectively. We observe an increase in the transmission as the field 
amplitude is increased due to intraband absorption bleaching starting at fields above 8 kV/cm. This 
effect arises from a field-induced reduction in THz conductivity that depends strongly on the Fermi 
energy. We account for intraband absorption using a free carrier Drude model that includes neutral 
and charged impurity scattering as well as optical phonon scattering. We find that although the 
Fermi-level dependence in the monolayer and five-layer samples is quite different, both exhibit a 
strong dependence on the field amplitude that cannot be explained on the basis of an increase in 
the lattice temperature alone. Our results provide a deeper understanding of transport in graphene 
devices operating at THz frequencies and in modest kV/cm field strengths where nonlinearities 
exist. 
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I. INTRODUCTION 

The well known linear electronic dispersion of graphene 
naturally gives rise to unusual optoelectronic properties,® 
potentially useful for high-speed modulators and the 
other active components. While it has been touted as 
an important material for terahertz (THz) frequency 
devices,® charge transport properties at THz frequencies 
under field conditions commonly found in devices, on the 
order lO’s kV/cm, are poorly understood. This is in part 
due to a lack of control over the extrinsic properties of the 
devices and graphene’s sensitivity to its environment.® 
The electrodynamic response of doped graphene has pre¬ 
viously been studied in the low field, linear regime under 
gated conditionJ^ for frequencies ranging from near dcP 
to THz,®IISl mid irP and optical regimes.^ Under in¬ 
tense fields, however, charge carriers can acquire signifi¬ 
cant momentum and for Fermi energies {Ef) close to the 
charge neutrality point (CNF), a complicated interplay 
between intra- and inter-band transitions can occurA 
better understanding of these nonlinear interactions and 
their density dependence is needed for the future device 
design, which is the aim of this study. For this purpose, 
we study the THz electrodynamic response of graphene 
under high field conditions and controlled Fermi level 
through a gate voltage Vg applied to an ionic gel top 


gate which is transmissive to THz light (see Fig. [^. 

The development of efficient nonlinear optical methods 
for the generation of intense THz pulses^ has made it 
possible to time resolve the high field response of charge 
and spin excitations in materials.!^ In graphene, due to a 
strong interband dipole matrix element that diverges at 
the charge neutrality (Dirac) point, the THz light-matter 
interaction can be easily pushed into the nonperturba- 
tive regime even by mode st electric fields on the order 
of several lO’s kV/ cm.E^^ As the thermalization time of 
accelerated hot carriers is noted to be exceedingly short 
in graphene,!!^ these nonlinearities can be quantified as 
an elevated temperature of the electronic subsystem.!^ 
Recent works have demonstrated a strong THz field- 
induced transparency in room-temperature graphen^^in 
both dopecP^HUl ^^d optically excited graphene.!^ This 
phenomenon has been described by ultrafast intra-band 
carrier thermalization on the order of 30 fs,^ leading 
to a transient elevated carrier temperature on the order 
of 2000 K and subsequent optical phonon emission dur¬ 
ing cooling. In this model, the non-equilibrium heating 
and generation of optical phonons subsequently enhances 
the carrier scattering, suppressing the field-driven current 
and thus THz absorption in graphene. Other sources of 
carrier scattering that can potentially play an important 
role in the THz response of graphene include long-range 
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FIG. 1. (a) Schematic of the graphene sample on a substrate 
with gold four-point-probe contacts and a side gate contact. 
The broadband THz pulses traverse the sample at normal 
incidence, (b) The mechanism for ionic gating of graphene is 
shown in the inset, illustrating the accumulated Debye layer 
under a gate voltage Vg- 


impurity scattering, due to the presence of the charged 
impurities in the substrate, short-range scattering, due 
to neutral impurities or disorder in the graphene itself, 
and carrier-carrier scattering. While impurity scattering 
dominates at low lattice temperatures, optical phonon 
scattering becomes more dominant at higher tempera- 
tures (T > 700 K)P 


Saturable absorption due to interband contributions 
has also been observed in the visible regime due to the 
phase-space filling effects as well as carrier heating. At 
THz frequencies, the threshold for an interband response 
is determined by the carrier density, as the carriers must 
be excited at twice the quasi-Fermi energy to undergo 
a transition. In this work, we study the Fermi energy 
dependence of the nonlinear THz transmission of gated 
graphene and show that although the nonlinear response 
is similar in different samples with different doping lev¬ 
els, it does not appear that this effect can be described 
by heating of the lattice alone. The paper is organized 
as follows. In Sec. H, we discuss the possible scattering 
mechanisms and their effects on THz transmission as a 
basis for the model used to describe the data. The results 
for both five-layer and single-layer graphene samples are 
presented in Sec. HI. Finally, the conclusions are sum¬ 
marized in Sec. IV. 


II. THEORETICAL MODEL 


In this section, we discuss the model that we use to 
account for the effects of tuning the Fermi energy on 
the transmission of THz light. Using this model, we 
extract the dependence of the scattering time on the 
Fermi energy and the lattice and electron temperature. 
There are four potentially important sources of carrier 
scattering in graphene: short-range impurity scattering, 
long-range impurity scattering, phonon scattering, and 
carrier-carrier scattering. 

Strictly speaking the carrier dynamics should be mod¬ 
eled using a full dynamic model that includes the nonlin¬ 
ear response arising from the linear band dispersion and 
inter band transit ions We find however, in agreement 
with our previous work,l^ that for the samples and field 
strengths considered here, there is essentially no intrin¬ 
sically nonlinear response and a linear Drude model is 
sufficient. The nonlinear response of the graphene is in¬ 
corporated in the model by employing a scattering time 
that depends on the THz field amplitude. Thus, in what 
follows, we employ the Drude model. 

First, we need to relate the transmission level to the 
conductivity of graphene. It can be easily shown that the 
transmission through the graphene sample normalized to 
the transmitted field through the substrate is given 


T( i 

Ts ~ l + n + Zoa{u)N' 


( 1 ) 


where Tq+s and Ts are the transmission through the 
sample and the substrate only, respectively, n is the re¬ 
fractive index of the substrate, aiuj) is the complex ac 
sheet conductivity of a single layer of graphene and N is 
the number of graphene layers in the sample. The con¬ 
ductivity, dominated by the intra-band response at THz 
frequencies is given by a simple Drude model 


(t(w) = 


D 


(1/r — iuj) ’ 


( 2 ) 


where uj is the angular frequency, r is the total carrier 
moment um scattering time, and D is the Drude weight 
given b}®^ 


D = 


2e^kBTe 

(irh?) 


In 


2cosh 


2kBTe 


( 3 ) 


where e is the electron charge, kB is the Boltzmann 
constant, Tg is the electron temperature, h is the re¬ 
duced Planck’s constant, and /i/ is the chemical poten¬ 
tial. For simplicity, the chemical potential is considered 
to be equal to the Fermi energy at room temperature. 
However, as the electron temperature rises, we adjust 
the chemical potential so that the net charge on the 
graphene is unchanged as the temperature rises. The 
adjusted chemical potential along with the correspond¬ 
ing temperature are used then to calculate the Drude 
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weight. Using Eqs. ([T]|^, we can extract the total ef¬ 
fective scattering time r for any measured transmission 
level for each gate voltage and THz field amplitude. The 
gate voltage and field amplitude dependence of the scat¬ 
tering time can therefore be directly obtained from the 
THz transmission data once the electron temperature is 
determined. 

We now turn to the model that we employ for the 
scattering time used in the Drude model. We present 
a simple model of the effects of the lattice temperature 
and the carrier density on the scattering times based on 
three scattering mechanisms: short-range neutral impu¬ 
rity scattering, charged long-range impurity scattering 
and scattering due to absorption of optical phonons. Al¬ 
though there can be contributions to the scattering rate 
due to interactions with acoustic phonons, these are usu¬ 
ally considered to be relatively small at the temperatures 
considered here.^^We do not include carrier-carrier scat¬ 
tering due to the difficulties in accurately modelling its 
effects and because, for similar Fermi energies and THz 
field amplitudes, it has been shown that carrier-carrier 
scattering can be neglected without qualitatively chang¬ 
ing the results.l^ We will return to this issue when dis¬ 
cussing our field-dependent results. 

The expression for the total scattering rate is thus 
given by 


1 1 1 

- —- ^ -, 

'^imp '^op 


(4) 


where Timp is the impurity scattering time Top is the opti¬ 
cal phonon scattering time. The impurity scattering rate 
is given by 


1 _ 1 1 

'^imp '^sr '^Coul 


1 


(5) 


In this expression, l/rgr is the short-range neutral impu¬ 
rity scattering rate, which has been shown theoretically 
to be linearly dependent on the energy of the carriers,!^ 
in our model, we take this rate to be proportional to the 
average energy per carrier with respect to the Dirac point 
K, which is given by 

E {Pee(K + k) + phhi'K + k)} hVF^i 

^ E{pee{K + k) + Phh{K + k)} ’ 

k 


where pee{phh) is the free electron (hole) population den¬ 
sity, and r’i? is the Fermi velocity. The second term in the 
impurity scattering rate is, rcouh which is due to charged 
long-range impurities. This rate has been shown to be 
inversely proportional to the square root of the carrier 
density (due to carrier screening)!^. Finally, Tmi is an 
additional gate-voltage-independent scattering time that 
we find is necessary in order to model the bias depen¬ 
dence of the scattering time in our multi-layer sample. 
The physical origin of this term is not clear, but it may 
be due to scattering from the ionic gel or is perhaps arises 


due to a correction term in the relationship between the 
gate voltage and the Fermi energy. Thus, the final ex¬ 
pression for the impurity scattering rate is 

4 _ „ C’ , ^Coul , 1 

- — Csr-t^av ^ -^-5 Dj 

'^imp -^eff '^ml 


where Csr and ccoui are constants that are experimentally 
determined and are dependent on the impurity densities, 
and Eeff is the effective Fermi energy that corresponds 
to the square root of the carrier density and is given by 
DtT?/^. 

We now turn to the scattering due to optical phonons. 
There are two potential types of optical phonons that can 
contribute to scattering: optical phonons in the grap hene 
itself and polar surface phonons in the substrate.^^EH 
We find, using the parameters in Ref. [4T] that, due to 
the weak coupling strength, the scattering rate contribu¬ 
tion from the surface phonons in both substrate materials 
is much smaller than those due to the graphene optical 
phonons. Thus, in what follows, we neglect the phonons 
in the substrate and include only the contribution from 
the phonons in the graphene itself.^^ 

Previous studies have shown that energy transfer from 
electro ns to phonons takes place on a timescale of about 
100 fsPI^ Thus, in principle, if there are electrons 
that have been driven by the field more than an optical 
phonon energy above any vacant states, one would ex¬ 
pect emission of optical phonons over the pulse duration, 
which would increase the lattice temperature. However, 
without performing a detailed simulation that tracks the 
energies of individual electrons, it is very difficult to accu¬ 
rately include this contribution to the scattering. Thus, 
for simplicity, we only include the process of phonon ab¬ 
sorption, but not phonon emission. 

The optical phonon scattering rate due to optical 
phonon emission in the graphene is given by 

1 ^ -Pq _ {Egy + huJo) 

Top {2ffVpamU)o) (e^ol^BTi _ 1) ’ 


where fiuOo = 147 meV is the graphene optical phonon 
energy, = 7.6 x 10“^ g/cm^ is the 2D mass density 
of graphen^^, Dq = 5 x 10^ eV/cm is the deformation 
potential for the optical phonons^ and Ti is the lattice 
temperature. Various values for the optical phonon en¬ 
ergy and defor mation potential have been reported in the 
literatur j^^ l ^^^ . We employ those found in Ref. [24l al¬ 
though, as we discuss in the following sections, we have 
considered the effect of modifying the deformation po¬ 
tential on the results. We finally note that in all cases we 
assume equality of the carrier and lattice temperatures 
that appear respectively in the Drude weight (Eq. (§) 
and the optical phonon scattering rate (Eq. (§)■ 
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III. EXPERIMENTAL RESULTS AND 
DISCUSSION 

In order to investigate the effect of the substrates and 
number of graphene layers on the transport properties 
of graphene to the intense THz field, we fabricated two 
different samples treated in the next two sections. The 
first is a five layer graphene sample, fabricated by trans¬ 
ferring sequentially monolayer graphene onto a sapphire 
substrate, and the second a monolayer sample on a quartz 
substrate. The substrates were patterned with 5/100 nm 
thickness of Cr/Au (prior to transfer) for the electrical 
four-point-probe measurements as shown in Fig. The 
gate was created by placing the fifth contact at the edge 
of the substrate so that it is in contact with the ionic gel 
for electrochemical doping. The ionic gel was fabricated 
using the method proposed by Ref. [35] and transferred to 
the sample under an inert glovebox environment. A uni¬ 
form layer of ionic gel with a thickness of less than 50 /im 
was spin coated on the graphene sample as well as the 
bare substrate to serve as a reference. Upon applying the 
gate voltage, mobile ions are transferred to the surface of 
the graphene to form a Debye layer with a thickness on 
the order of 1 nm.l^This type of solid polymer electrolyte 
gate is much thinner than the conventional 300 nm Si02 
back gate, thus a much higher gate capacitance can be 
achieved. With the electrolyte ionic gating, high doping 
levels on the order of 5 x 10^^ cm“^ are possible.l^ The 
gating response of the graphene is measured in terms of 
the electrical current driven from the source to the drain 
electrode through the graphene channel as a function of 
Vg- In order to estimate the Fermi energy as a function of 
gate voltage relative to the charge neutrality point volt¬ 
age (Vgatp), we use the following equation from Ref. [101 
for a similar geometry and ionic gel, where 

Ef = 0.346y/ {Vg — VcNp) [^y]- (9) 

The case of Vg = Vcnp needs to be treated sep arately 
because, in agreement with other resear chers,[^ we find 
that due to spatial inhomogeneity of the charge density 
across the graphene sample an effective Ef must be de¬ 
termined at the nominal CNF gate voltage, in a range of 
±100 meVP 

The high-held time-domain THz spectrometer was de¬ 
veloped based on tilted-pulse front optical rectihcation in 
lithium niobate,^ with peak THz held amplitudes lim¬ 
ited to 70 kV/cm in this work. The bandwidth of the 
THz pulses was ^ 3 THz and was detected using electro¬ 
optic sampling in a 300 /im thick (llO)-cut GaP crystal. 
The emitted THz pulses were collected and focused by a 
set of oh-axis parabolic mirrors to a 450 /rm spot size as 
determined by an uncooled microbolometer camera. The 
graphene sample was positioned carefully at the THz fo¬ 
cus for high-held transmission measurements at normal 
incidence. The amplitude of the THz pulses was varied 
using a pair of wire-grid polarizers prior to the sample in 
the collimated section of the beam. The measurements 
were performed by sweeping the graphene gate voltage 
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FIG. 2. (a) Four-point-probe measurements of the 5-layer 

sample dc resistance with the charge neutrality point indi¬ 
cated by the resistance peak, (b) The Raman data (with 
the G and 2D peaks) of single layer of CVD grown graphene 
before stacking to form the 5 layers. 


and registering the value of the transmitted peak THz 
electric held. No THz pulse reshaping was observed as 
the gate voltage varies from —2.5 V (below the Dirac 
point) to 2.5 V, thus all information on the nonlinear 
transmission can be quantihed by monitoring the THz 
pulse peak held transmission. 


A. Five-layer sample 

Large-area, single-layer graphene was fabricated by 
chemical vapor deposition (GVD) on Cu-foil using stan¬ 
dard techniques. Poly(methyl methacrylate) (PMMA) is 
used to transfer the graphene to a c-cut sapphire sub¬ 
strate with refractive index of 3.31 and the PMMA layer 
was removed. The process was repeated hve times so that 
the hve-layer graphene sample was obtained, with the 
layers electrically decoupled from one another. Fig. da) 
shows the sheet resistance of graphene as a function of 
gate voltage. The GNP is specihed as the resistance 
maximum which occurs in this case at Vg = —0.8 V. 
We note that hysteresis was observed in these resistance 
measurements when sweeping the gat e volt age, which is 
commonly observed in ionic gel gates.l^^^^ We therefore 
only perform measurements on the positive sweep of the 
gate voltage for repeatability. The Debye absorption of 
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FIG. 3. For the five-layer sample: (a) The THz pulse peak 
held transmission through the gated graphene sample for the 
held strengths indicated, normalized to the transmission at 
the Dirac point at Vb — Vca^p = 0. (b) The held dependence 
of the THz peak held transmission | t |= various 

gate voltages. Einc is the incident THz electric held at the 
graphene sample position. 


the ionic gel is located at lower frequencies than THz, 
producing a transparent top gated graphene sample. All 
measurements in this work are performed at room tem¬ 
perature. Fig.j^b) shows the Raman spectroscopy of one 
of the hve layers of graphene before stacking. The ratio 
between the 2D and G peak values is 2.8 which indicates 
high quality single layer graphene .SSI 

Fig. [^a) shows the variation of the transmitted peak 
THz electric held (normalized to its maximum value at 
the Dirac point) as a function of gate voltage sweep 
for different incident peak electric helds in the range 
of 2 kV/cm to 70 kV/cm. The THz modulation in¬ 
duced by the gate voltage sweep (Fermi level change in 
graphene) is a strongly nonlinear function of the peak 
electric held value. At low helds, gate-induced THz mod¬ 
ulation is more than 31% whereas for higher held ampli¬ 
tudes the modulation decreases and saturates at a value 
of 7.4% above 60 kV/cm. As in previous work on un¬ 
gated graphene, a saturable power transmission function 
can be used to phenomenologically describe the behavior 
observed in Fig. [^b).II^ The power transmission is de- 

hned as | T |= ^ ^ ^ where Egamp is the peak THz 

electric held transmitted through the ionic gel, 5 layers 
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FIG. 4. For the hve layer sample: (a) THz power transmission 
dehned as | T |= curve htted by Eq. 10 for two 

^ref 


diherent gate voltages, one at the Dirac point (Vg — Vcnp = 
0 y, red line) and the other one at a n-doped Fermi level 
corresponding to (Vg — Vgnp = —2.4 E, blue line), (b-d) Fit 
parameters Esat, Tun, and Tns extracted from Eq. (10) for 


power transmission curves as a function of gate voltage. 


of graphene and the substrate, and E^ef is the peak held 
transmitted through the ionic gel and the substrate only. 
The model is given byf^^ 


T{Einc) = T; 


’ TTi 2 / ^ 2 

J^inc / J^sat 


( 10 ) 


where Ei^c^ Egat^ Tn^ and T^g are the incident peak elec¬ 
tric held, the saturation electric held, the linear power 
transmission coefficient and the nonsaturable power 
transmission coefficient, respectively. 

In Fig. [^a) we present the measured transmitted 
power as a function of the incident held along with the 
ht using Eq. ( [Tq| ) for Vg = Vcnp (corresponding to 
the CNP, red line) and Vg — Vcnp = —2.4 V (i.e. 
highly the n-doped case, blue line). The saturation on¬ 
set of the power transmission, |T|, occurs at a lower 
Vine = 9.7 kV/cm at the CNP relative to what occurs 
for higher doping, as demonstrated by Fig. [^b). Tun 
[Fig. gc)] and Tns [Fig. |4];d)] are higher at the CNP 
than at other Fermi levels, overall expected due to the 
larger conductivity of the doped graphene. The satura¬ 
tion behavior for the curves in Fig.|^a) is explained by a 
decrease in carrier mobility as the THz excitation redis¬ 
tributes carriers within the conduction band.^ There¬ 
fore, the decrease in Egat when the gate voltage ap¬ 
proaches the Dirac point can be attributed to an increase 
in the scattering rate when the Fermi level is tuned to the 
CNP, as we shall discuss shortly. 

We now use the model presented in Section II to ex¬ 
tract the scattering time as a function of the Fermi energy 
and the incident THz field. While the Fermi energy may 
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not be precisely the same in all 5 layers of the graphene 
sample, for simplicity we use the same value for the five 
layers. Similarly, we assume that the scattering time is 
the same in all layers. We estimate Ef from gate voltages 
other than the CNP voltage using Eq. To determine 
the Fermi energy for Vq = Vcatp, we first use the model 
of Section II to extract the fitting parameters, ccoui 
and Tmi in Eq. for the transmission data for the lowest 
THz field amplitude of 2 kV/cm. Taking the lattice and 
electron temperatures at this lowest field to be 300 K, we 
obtain Cs^=49 ps“^eE“^, ccoui = 0 and =7.3 ps“^. 
The data as well as the fit are shown in Fig. j^a). We now 
extract an effective Fermi energy at the Vg = Vcnp of 
approximately 90 meV by extrapolating the fitting curve. 
This value is in good agreement with other studies that 
used the same mechanism to tune the Fermi energy.lISl We 
find from our fit that the contribution of charged impu¬ 
rities is negligible and that at this temperature and field 
strength, the contribution of the optical phonons is also 
very small. Thus, the dominant contribution to scatter¬ 
ing at this lowest field appears to arise from short-range 
impurities, which results in the approximately hyperbolic 
relationship between r and Ef seen in the figure 

We now turn to the dependence of the scattering time 
on the incident field for different Fermi energies. Rather 
than using our temperature-dependent model of scatter¬ 
ing, we first assume that the electron and phonon tem¬ 
peratures do not change with the incident THz field am¬ 
plitude and simply extract the net scattering time as a 
function of the incident field strength to two different 
Fermi levels. This is plotted in Fig. ib). As can be 
seen, for both Fermi energies, the scattering time de¬ 
creases dramatically with increasing field amplitude. We 
see that the decrease in the scattering time with field is 
much larger for the low Fermi energy case. This would 
be expected if the increased scattering is due to a mech¬ 
anism that has a rate that is roughly proportional to 
the Fermi energy and that would increase at high fields. 
There are two likely candidates: optical phonon scatter¬ 
ing and carrier-carrier scattering. 

We first consider optical phonon scattering. When the 
field is higher, it heats the graphene up considerably; this, 
in turn results in a larger optical phonon population and 
a shorter scattering time. To test the viability of the 
optical phonon mechanism, we employed our model that 
includes the effects of electron and lattice temperature. 
Using the impurity scattering rates found at the lowest 
field and requiring that the electron and phonon temper¬ 
atures are equal, we were able to obtain good agreement 
with the experimental data for the transmission by using 
the temperature as a fitting parameter. We found that at 
the highest fields, the temperature required was approx¬ 
imately 1550 K. To determine if such temperatures are 
physically sensible, we compared the experimental result 
for the energy absorbed per unit area to the increase in 
the electron and phonon populations per unit area. We 
found that the required increase in the system energy 
was more than an order of magnitude greater than the 


absorbed THz energy for the 70 kV/cm field. Most of 
the increase in the energy of the system when the tem¬ 
perature increases to such levels resides in the optical 
phonons, not the electrons. Because there is not uni¬ 
versal agreement as the values of the optical phonon en¬ 
ergy and deformation potential, we tried adjusting these 
parameters. However, we could not find any values of 
these parameters that yielded anything approaching en¬ 
ergy balance at all fields and Fermi energies for one set of 
parameters. We conclude from this, that for this sample 
at least, the decrease in the scattering time with increas¬ 
ing field does not arise entirely from an increase in the 
optical phonon population due to an increase in the lat¬ 
tice temperature. 

Another possible mechanism for the field dependent 
scattering is carrier-carrier scattering. This scattering 
rate increase rapidly as the carrier density increases (as 
observed). Moreover, it also increases with the energy of 
the carriers, which will occur when the field amplitude 
is high. Modeling carrier-carrier scattering with energy 
dependence would require a sophisticated Monte-Carlo 
or density matrix simulation. Thus, although we cannot 
rule out this mechanism, it is beyond the scope of the 
current work. 

We now examine the results for our monolayer 
graphene sample to see how its transmission depends on 
the Fermi energy and the THz field amplitude. 


B. Single-layer sample 

A single-layer graphene sample was prepared by chem¬ 
ical vapor deposition on Cu-foil and transferred onto a 
z-cut quartz substrate with a refractive index n = 1.96. 
The sheet resistance of this sample as a function of Vg 
was determined by a two terminal measurement shown in 
Fig. I^a) with the CNP indicated again by the resistance 
peak. Electron (hole) doping is induced by changing Vg 
above (below) the CNP voltage Vgnp and Eq. |^is used 
to calculate the Fermi level energy. Fig. [^b) shows the 
Raman spectra for the single layer graphene sample, with 
a 2D/G peak ratio of 1.78 confirming high quality mono- 
layer graphene. 

Fig-iii: a) shows the normalized THz peak electric field 
amplitude transmitted through the graphene when the 
Fermi-level is swept from -1 to 0.5 V. Again, our data 
does not show any phase change in the THz waveform 
transmitted through the graphene by increasing the in¬ 
cident THz field and so the conductivity change can be 
completely quantified by monitoring the peak transmis¬ 
sion. For the lowest THz peak electric field strength of 1.5 
kV/cm, a significant 20% reduction in the peak transmis¬ 
sion was observed at Vg - Vgnp = —lVorEj = —340 
meV. As expected, the increased doping leads to higher 
conductivity and therefore absorption of the THz pulse, 
however, increasing the incident THz field leads to a re¬ 
duction in t he do ping-induced absorption as in the 5- 
layer samplc.^^^^^ As in the 5-layer case, a continuous 
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FIG. 5. (a) The impurity scattering time (r^mp) versus Fermi 
energy, and (b) total scattering time versus incident field am¬ 
plitude at the CNF and at the highest Fermi level of 442 meV 
for the five-layer sample. Note that in (a), the value for Ef 
the first point is determined via a fit to the curve. 


reduction in field-induced transmission modulation oc¬ 
curs as the gate voltage approaches the charge neutrality 
point. This behavior is best seen in Fig. [^b), showing 
the peak transmission as a function of the THz peak field 
strength for various applied Vg- The THz field-induced 
transparency (the percentage increase in transmission) 
is more pronounced for high doping levels. As with the 
five-layer sample, we find that the induced nonlinearity is 
mainly due to an increase in th e carrier scattering when 
the THz field is increased. 

We extract the scattering behavior by using the model 
of Section IT We extract the low-field scattering parame¬ 
ters in the same way as was done in the last section taking 
the temperature to be room temperature. The transmis¬ 
sion data at the lowest field amplitude can be well repro¬ 
duced by the theory only if we consider both short-range 
neutral and long-range charged impurity scattering in ad¬ 
dition to optical phonon scattering. We find that the best 
fit is obtained for Cgr ccoui =3.63p5“^eH 

and and Tmi =0. The data as well as the fit are shown in 
Fig. [^a). Note that the increase in the scattering time 
with Fermi energy indicates that long range charged im¬ 



FIG. 6. (a) Two terminal measurements of the resistance 

and (b) the Raman data (with the G and 2D peaks) of the 
monolayer sample. 


purity scattering is dominant at this field. The presence 
of long range scattering in this sample may be linked 
with the direct contact of the graphene with the Si 02 
substrate as opposed to the 5 layer sample, where 4 of 
the layers are physically removed from the substrate and 
the effects of charged impurities on carriers in the up¬ 
per layers are partially screened by carriers in the lower 
ones. Finally, using the extrapolation of the fit to low 
Fermi energy, we extract and Fermi energy of 88 meV 
for the bias Vg = Vgatp- 

We now use the same constant-temperature model of 
the previous section to determine the net scattering time 
as a function of the incident field strength for the five dif¬ 
ferent Fermi levels and plot these times in Fig. |^b). As 
in the five-layer case, the scattering time decreases dra¬ 
matically with increasing field amplitude for all Fermi en¬ 
ergies. However, in contrast to the results of the previous 
section, we find that the decrease in the scattering time 
with the field is much smaller when Vg = Vgatp- One 
reason for this difference is that in the monolayer, the 
low-field scattering time is the shortest when the Fermi 
energy is the lowest, thus the effect of any additional scat¬ 
tering channels is smaller. At the highest fields, just as 
in the five-layer sample, the scattering time is the largest 
when Vg = Vgatp- This seems to indicate again that 
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FIG. 7. (a) The normalized peak transmission as a function of 
the gate voltage for various THz peak electric field strengths 
and (b) the dependence of the peak transmission on the THz 
peak electric field for various doping levels for the monolayer 
sample. 


FIG. 8. (a) The impurity scattering time (rimp) versus Fermi 
energy at the lowest incident THz field of 1.5 kV/cm, and (b) 
the THz field dependence of (a) the carrier scattering rate. 
Note that in (a), the value for Ef the first point is determined 
via a fit to the curve. 


the field-dependent scattering rate is larger when the car¬ 
rier density (or perhaps average carrier energy) is larger. 
However, the trend is not the same for the intermediate 
Fermi energies, where at the highest field, the scattering 
rate is larger when the Fermi energy is smaller (with the 
exception of the CNF result). It appears that the result 
at the CNF is somewhat of an anomaly in this sample. 

Although, the Fermi energy dependence of the field- 
dependent scattering component is not as clear as in 
the five-layer sample, we again examined where optical- 
phonon scattering is a possible mechanism. We modeled 
the optical phonon scattering by taking the temperature 
to be a fitting parameter. Again, although we could re¬ 
produce the observed transmission curves, the tempera¬ 
tures required were very large. The temperatures at the 
highest field ranged from 1300 K for the highest Fermi en¬ 
ergy to almost 3000 K at the CNF. Again, such huge tem¬ 
peratures require that the energy deposited in the system 
is a factor of 3 to 10 larger than the energy that is experi¬ 
mentally found to be absorbed by the graphene. As in the 
five-layer case, we tried adjusting the phonon energy and 
deformation potential, but for all combinations that we 


tried we could not achieve anything approaching energy 
balance at all fields and Fermi energies. We thus conclude 
that for this sample as well, optical phonons cannot be 
the major contributor to the increase in the scattering 
rate with increasing field. In this sample, however, due 
to the more complex dependence of the field-dependent 
scattering contribution on the Fermi energy, it is not as 
clear that carrier-carrier scattering is a likely candidate. 

IV. CONCLUSION 

We have studied nonlinear THz light-matter interac¬ 
tion for both monolayer graphene on quartz and a stack 
of 5 monolayers of graphene on sapphire under the con¬ 
trolled Fermi level conditions. The two samples exhibited 
opposite dependencies on the Fermi energy for low THz 
fields. The monolayer sample scattering rate was domi¬ 
nated by charged impurities, while the five layer sample 
scattering was neutral impurity dominated. This differ¬ 
ence is likely explained as arising from different prepara¬ 
tion recipes as well as differences in the substrates. 

We also observed a strong increase in the THz trans- 
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mission with increasing field that occurs for THz peak 
electric fields as low as 8 kV/cm. This increase appears 
to result from a decrease in the carrier scattering times 
with increasing field. Our analyses indicate that opti¬ 
cal phonon absorption is not likely the main source of 
this field-induced increase in the scattering rate. The 
source of this increase remains uncertain. For the five 
layer sample, the most promising mechanism is perhaps 
carrier-carrier scattering, but for the monolayer sample 
the mechanism is even more unclear. 

This work points to the importance of future modelling 
of full carrier and lattice dynamics in graphene excited by 
intense THz pulses to determine the main source of the 
strongly nonlinear response. The interaction of high field 


THz pulses with graphene is significantly dependent on 
Fermi energy. Moreover, given the variability in process¬ 
ing of graphene films, we show that it is entirely possible 
that different impurity scattering mechanisms dominate 
in different samples, even though the field dependence is 
similar for a fixed Fermi energy. Thus any comprehen¬ 
sive study must be under controlled gated conditions. We 
anticipate that our findings will be important for future 
designs of graphene devices operating at THz frequencies, 
influencing performance due to nonlinearities arising at 
the kV/cm operating fields. 

The authors gratefully acknowledge funding support 
from NSERC, FRQNT and CFL 
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